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Abstract. 

The spectrum of glueballs in 2+1 dimensions is calculated within an extended class 
of Isgur-Paton flux tube models and compared to lattice calculations of the low-lying 
SU(iV > 2) glueball mass spectrum. Our modifications of the model include a string 
curvature term and a new way of dealing with the short-distance cut-off. We find that 
the generic model is remarkably successful at reproducing the positive charge conjuga- 
tion, C = +, sector of the spectrum. The only large (and robust) discrepancy involves 
the _+ state, raising the interesting possibility that the lattice spin identification is 
mistaken and that this state is in fact 4r + . Additionally, the Isgur-Paton model does 
not incorporate any mechanism for splitting C = — from C = + (in contrast to the 
case in 3+1 dimensions), while the 'observed' spectrum does show a substantial split- 
ting. We explore several modifications of the model in an attempt to incorporate this 
physics in a natural way. At the qualitative level we find that this constrains our choice 
to the picture in which the C — ± splitting is driven by mixing with new states built 
on closed loops of adjoint flux. However a detailed numerical comparison suggests that 
a model incorporating an additional direct mixing between loops of opposite orienta- 
tion is likely to work better; and that, in any case, a non-zero curvature term will be 
required. We also point out that a characteristic of any string model of glueballs is 
that the SU(iV — > oo) mass spectrum will consist of multiple towers of states that are 
scaled up copies of each other. To test this will require a lattice mass spectrum that 
extends to somewhat larger masses than currently available. 



PACS numbers: 11.15.Kc, 12.38. Gc, 12.39.Mk 



1 Introduction 



While it is now possible to calculate the spectrum of continuum non-Abelian gauge 
theories with some precision, using standard lattice Monte Carlo techniques []!], |2|, [|, 
we know little about the structure of these glueballs. This is to be contrasted with 
states containing quarks where, at least for the low-lying spectrum, the quark model 
provides a remarkably successful semi- quantitative model framework for understanding 
the structure of mesons and baryons. (Apart from the interesting cases of scalar mesons 
and pseudoscalar flavour-singlet mesons.) 

The glueballs of the SU(3) non-Abelian gauge theory in 3+1 dimensions are par- 
ticularly important because the presence of such extra non-quarkonium states in the 
spectrum of QCD (and in the experimental spectrum) would provide a direct reflection 
of the gauge fields in the theory. Understanding just how they do so (mixing, decays 
etc.) would be made easier if we understood something about their structure. Unfor- 
tunately, beyond providing some information about glueball sizes, lattice Monte Carlo 
calculations have as yet given us little insight into their structure. Such calculations 
involve connected correlators of several operators and while they are simple in princi- 
ple it is, in practice, much harder to achieve sufficient statistical accuracy than in the 
corresponding mass calculations. 

An alternative way to learn about the structure of glueballs is through a reliable 
model - just as the quark model provides us with useful information on the structure 
of the low-lying mesons and baryons. To establish whether a glueball model is 'reliable' 
one can compare the spectrum it predicts to the known spectrum (as calculated from 
the lattice). This is the approach we follow here. There are two obvious models that 
one might try: constituent gluon models, such as gluon potential || and bag models 
0, or flux tube (string) models 0. In this paper we shall confine ourselves to a study 
of the latter. 

Models try to isolate the essential physics and neglect everything else; thus they 
necessarily involve approximations. So one does not expect precise agreement with 
the known spectrum. If we are only looking for semi-quantitative or even qualitative 
agreement, it is important to test the model in as many relevant contexts as possible. 
One fact we can usefully use here is that string models (and indeed bag models) can be 
equally motivated in any gauge theory that has linear confinement and hence string-like 
flux tubes. This suggests that it would be useful to test this model not only in 3+1 
dimensions |7j but also in 2+1 dimensions where non-Abelian gauge theories appear 
to be linearly confining and detailed mass spectra are available ||. This is what we 
shall do in this paper. (We remark that recent improvements in the lattice calculation 
of the D=3+l SU(3) spectrum warrant a complete update of the study in 0.) 

In the next section we briefly review the Isgur-Paton flux tube model for glueballs || 
and describe the qualitative features of the mass spectrum that it predicts for SU(N) 
gauge theories in D=2+l. We compare this spectrum with the 'true' spectrum as 
calculated on the lattice, |J, which, for the reader's convenience, we summarise in 
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Table [l] borrowed from |§ . We point out where the main discrepancies and difficulties 
lie and we point out some compelling generalisations of the model. In the subsequent 
section we address a major such difficulty: how to incorporate an acceptable C — ± 
splitting into the model. We finish with a summary of our results. Finally we remark 
that a summary of some of our preliminary results has appeared elsewhere f9j, and a 
much more detailed exposition appears in Jl0| . 



2 The Isgur-Paton flux tube model of glueballs 

Consider a quark and an antiquark sufficiently far apart. In a linearly confining theory, 
they will be joined by a flux tube which contributes an energy that is approximately 
proportional to its length. One can attempt to use such states as the basis for a flux 
tube model of quarkonia. The corresponding model for glueballs would be based on a 
loop of fundamental colour flux that closes on itself. This colour singlet object contains 
no quarks. If we neglect its thickness, we have a closed string of flux, and the glueball 
mass spectrum is obtained by finding the energy eigenstates of the quantised string. 
This is the starting point for the Isgur-Paton model M and our consequent extensions. 

We start with a closed string of colour flux in the form of a circle of radius p with 
bare string tension a b , and hence a bare energy 

E b = 2na b p. (1) 

Fluctuations about this circle are decomposed into discrete phonons of definite helicity. 
These phonons carry angular momentum, and so there must be a term proportional to 
the total phonon number M, given by 

- M - 1 ST ( + 

-^phonons = / . ^H^m ^m) \"J 

P 9^2 

added to the total energy of the state. Note that the above sums begin with the m = 2 
mode. The m = 1 mode is excluded in the model || because infinitesimal m = 1 
fluctuations are the same as infinitesimal translations of the circle. 

As usual when one quantises over modes of all frequencies there is a divergent con- 
tribution to the vacuum energy. Part of this divergent piece can be absorbed into the 
renormalised string tension, o = o b + c/2tt. The rest appears as a string Casimir en- 
ergy. It is universal for bosonic strings and, in the case where the strings end on static 



quarks, is called the Luscher term |LT|. In the present case the string has periodic 
rather than fixed boundary conditions, just as in calculations involving Polyakov loops 
[Oj, but the coefficient differs due to the exclusion of the m = 1 mode. Putting all 



this together we can write the energy of the string plus its modes as 

13 1 

E s = 2npo - — + - m ( n m + U m) ( 3 ) 

LZ P P m=2 
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However, we know that what we have is not really a one-dimensional string but rather 
a flux tube whose width will be ~ l/V "- For the low- lying states of interest to us we 
expect p ~ 1/ -y/er and so one might expect the simple harmonic fluctuations of the flux 
tube to be somehow suppressed. This is incorporated within the Isgur-Paton model 
|| by multiplying the contribution of the phonons to the energy (including that of the 
vacuum) by a heuristic suppression factor that rapidly approaches unity as p increases. 
This leads to a final string energy 

E?(p) = 2*pa + ¥^lF(p) (4) 

where 7 = —13/12, M is defined in eqn(g) and F(p) is the factor that suppresses 
the string excitations at small p. In the original model this was chosen to be F(p) = 
1 — e~^ p , where / is a parameter which we would expect to be 0(*Ja). This form is 
reasonable but somewhat arbitrary; one might ask, for example, why the string energy, 
27rp<r, is not modified at small p as well. 

To quantize the string, we must take into account its motion in the radial direction 
as well. In the Isgur-Paton model one identifies the conjugate momentum for the string 
and writes a Schrodinger equation in the radial coordinate p: 

where £ = p 3//2 turns out to be the natural variable to use here. This formalism 
assumes that the phonon modes are 'fast' compared to the collective radial modes and 
that they can therefore be treated as providing an effective potential for these latter 
modes. Clearly such an 'adiabatic' assumption is at best approximate: the model has 
only one scale, y/a, and so there is no reason for the phonon fluctuations to be much 
faster than the collective radial fluctuations for the low-lying part of the spectrum that 
will interest us. (Indeed if one calculates the low-lying spectrum one finds that the 
energy splitting associated with an increment in the phonon number is of the same 
order as the splitting associated with an increment in the radial quantum number. 
This suggests that the division into fast and slow modes is a crude approximation at 
best.) 

If we were in 3+1 rather than in 2+1 dimensions, the above description of the model 
would change as follows. First, rotations of the flux-loop around a diameter provide 
an additional source of angular momentum, and equation^) acquires a corresponding 
angular momentum term. In addition there are extra phonons arising from fluctuations 
of the loop normal to its plane. This doubling of modes leads to a doubling of the value 
of the string Casimir energy in eqn(|3|). 

The simplest version of the Isgur-Paton model sets F(p) = 1 in equation^) so that 
there is no fudge factor. Since a merely sets the overall scale of the mass spectrum, the 
mass ratios mj sfd are then predicted with no free parameters at all. These prediction 
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(borrowed from reference 0) are listed in Table [|. The comparison is with SU(N — > 
oo), since in that limit, just as in the model, heavy glueballs do not decay. The overall 
qualitative agreement is remarkable, with only the 0~ + far from its prediction in the 
C = + sector, and motivates the more detailed investigation of this paper. 



3 Generalising the Isgur-Paton flux tube model 

In this section we point to several ways in which the original Isgur-Paton model can 
be generalised. We start with the observation that one can build on other strings 
than the fundamental. We then point out that a curvature term in the effective string 
action makes an important difference. The argument for both these extensions is 
compelling. We then turn to the short-distance fudge factor F(p) in eqn([|), point 
out its shortcomings and suggest some alternatives. We leave to the next section the 
important question of how to split the C = ± spectra in a way that is both natural in 
terms of the string model and reproduces the main features of the observed splitting. 



3.1 Extra strings, extra states 

The flux tube in the Isgur-Paton model contains flux in the fundamental representation; 
it joins charges that are in that representation. For SU(iV > 4) there exist charges in 
higher representations which cannot be screened by virtual adjoint charges (i.e. gluons) 
down to the fundamental. One can label charges in these representations by the way 
they transform under a centre gauge transformation, z G Z^. If they acquire a factor 
z k we will refer to them as having iV-ality k. Since gluons transform trivially under 
the centre they cannot screen the iV-ality of a charge. For each such charge we have 
a flux tube of a corresponding iV-ality k, which will possess a string tension <jfc. We 
can consider a closed tube of such flux, and we can then build a whole spectrum of 
glueball states on this flux string just as we did for the fundamental, k — 1, string 
in the Isgur-Paton model. Thus as iV grows the spectrum will acquire extra towers 
of states that are identical to the spectrum obtained with the fundamental flux loop 
except that their overall energy scale is \jo^Jo\ (ignoring any mixing). 

If observed, such a spectrum would be a remarkable manifestation of the underlying 
string structure of glueballs. Of course it is not guaranteed that such a spectrum 
actually exists in the string picture: this will depend on the dynamics. Consider for 
example the case of SU(4). If it happens to be the case that Ok=2^^k=\ then the 
k = 2 flux tube can break up into two k = 1 (fundamental) flux loops, so that the 
k = 2 states are just multi-glueball scattering states formed out of the k — 1 glueball 



states. In the SU(4) gauge theory it is known jnj that a k=2 ~ 1.4cr fc=1 in both D=2+l 
and D=3+l, and so the corresponding extra states should exist there. However the 
lattice calculations have not identified enough excitations (in each J PC sector) to test 
for the possible presence of such extra states, so we will ignore this potential state 
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replication in the remainder of this paper, apart from pointing to its great interest for 
future calculations flQfl . 



3.2 Curvature/elasticity 

The flux tube must have a finite thickness if it is to have a finite energy density; 
presumably it will be 0(1/^0"). Such a finite flux tube will presumably possess an 
effective elasticity. In string language this is a curvature term. For a meson, the 
curvature of the straight string joining the quarks is zero and so a curvature term 
would have no effect to the order in 1/p that we are including. For a closed string, on 
the other hand, the curvature is constant and integrates to a ~ 1/p contribution. The 
constant of proportionality, our effective elasticity, we denote by 7^ and we will regard 
it as an unknown free parameter. Note that we may regard the Casimir energy of the 
closed loop as simply renormalizing the elasticity 

1 = 1e-— (6) 

just as the cp piece was absorbed into a renormalisation of o. Although there has been 
some discussion as to the sign such an elasticity should take, we shall leave 7 as a 
free parameter, whose value is to be determined by fitting the spectrum. 



3.3 Modification at short distances 

Since the flux tube has a finite width, a glueball will presumably cease to look like an 
excited closed string when p is much less than that width. This is embodied in the Isgur- 
Paton model by a fudge-factor F(p) = l — e~f p which suppresses the contribution of the 
string phonons as p — > 0, as in eqn(^). The detailed form of F(p) is largely arbitrary, 
as is the choice to suppress the phonon excitations but not the 2np string contribution. 
Since the spectrum of the string model is non-singular when we set F(p) = 1, the effects 
of the suppression factor are not large and the details do not matter greatly. We have 
calculated the spectrum for various possibilities and we find that as far as the C = + 
spectrum is concerned what one needs is a modest short-distance suppression so as to 
get the ++ ,2 ++ splitting about right, and then one can tune the 7/p contribution so 
as to raise the overall spectrum to about the right level. 

A quite different possibility is to make the string tension a function of p rather than 
to impose a fudge-factor F(p). This is motivated by a recent study |Tj| of closed flux 



tubes in the dual Ginzburg-Landau theory. They find that the effective string tension, 
a e ff(p), varies with p so as to vanish as p — > 0. One can in fact parameterise the 
observed variation quite accurately using 



~{p) = a{l-e-^). (7) 
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One can then quantise the string model with this p-dependent string tension and solve 
for its spectrum. Since a e ff(p) appears in the mass and hence in the kinetic energy of 
the loop, the quantisation is not entirely straightforward and we leave its description 
to the Appendix. (More details may be found in ]n|.) The qualitative effect of using 
a variable a e jj(p) in eqn@ is that at small p the kinetic energy is enhanced relative 
to the potential energy. This is much the same as the effect of our fudge-factor F(p). 
However it has the advantage that there is no free parameter (or functional form) to 
choose and there is no ambiguity as to how one should apply it. 



4 Splitting C = + from C = — 

For N > 2 the flux tube carries an arrow. In the simple Isgur-Paton model there is no 
mixing between loops of opposite flux, and the resulting states are degenerate. Since 
the direction of the flux reverses under charge conjugation C, this predicts degenerate 
C = ± spectra, as in Table 0. We now turn to the problem of how one might split 
this C — ± degeneracy in a way that is both natural in terms of the string model and 
reproduces the main features of the observed splitting. We shall begin by summarising 
what these features are and we shall then consider two possible dynamical mechanisms, 
direct mixing and adjoint mixing. Additional mechanisms, such as indirect mixing and 
k-string mixing, are explored in flQ |. In each case we shall ask how well the main 
observed features are reproduced. 

4.1 The observed C = +/— splitting 

From the masses listed in Table [TJ we see that the C = +/— splitting possesses the 
following qualitative features: 

• In Fig.|l] we plot two examples of the C — ± splitting, as a function of 1/N 2 . This is 
a natural variable to use since the leading corrections to the large- N limit are expected 
to be 0(1/N 2 ) [|l3j. We infer from this plot (and from similar plots of other states) 



that the splitting remains non-zero in the N = oo limit; it is a leading order effect. 
• In SU(2) there is no C = — sector and one can ask what the SU(2) + mass continues 
to as iV increases. This will clearly depend on the dynamics that produces the C = ± 
splitting for iV > 3. For example, if this dynamics simply splits the ++ and 
equally from their naive degenerate masses, then we would expect the SU(2) + to 
continue smoothly to the average of the ++ and masses. As another example, if 
the shift involves just the , then the SU(2) + will continue smoothly to the N > 3 
++ . Conversely, if the shift affects just the ++ , then the continuation should be with 
the . In Fig.|| we plot the ++ and masses, as well as the average of the two, 
as a function of 1/N 2 . We see that in all cases the variation with N for N > 3 can 
be described using just a leading oc 1/N 2 correction. We also note that the SU(2) + 
mass extrapolates precisely from the ++ masses while it is inconsistent with a smooth 
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extrapolation of the averaged C — ± states or of the . The same is true for the 
tensor: the SU(2) 2 + mass is a smooth continuation of the SU(iV > 3) 2 ±+ masses, 
and not of the average of the 2 ±+ and 2 ±_ masses, or of the 2 . We shall see that 
this observation provides a tough constraint on possible mechanisms for splitting the 
C = + and C = — sectors. 

• The C = ± splitting appears to decrease as the mass increases. To be more specific 
we infer from 

m — — m ++ > m — * — m ++*. (8) 

that states with larger radial quantum number, tir, are split less than those with 
smaller hr. (Recall that in the flux tube model these lowest excitations are radial 
rather than phonon.) Furthermore we infer from 

m — — m ++ = 1.85(26)^/(7 > m 2 — — m 2 ++ = 1.01(38)a/<t 

> ma— - m 1+ + = -0.62(65)^/0- (9) 

(obtained at iV = oo) that the magnitude of the splitting also decreases with increasing 
phonon number. (Recall that for the lightest J = 0, 2, 1 glueball states the total phonon 
number in eqn(Q) is M = 0,2,5 respectively.) Note that the decrease we see is even 
faster when the splitting is expressed in terms of the average mass. All this provides 
constraints on possible splitting mechanisms. 



4.2 Direct mixing 



Since a flux loop has a direction, L or R, it is convenient to introduce 2-component 
wave-functions: 



% = 



(10) 



In this notation, we can write (unnormalized) C — ± states as 
and eqn([|) becomes 







H IP V 













Hr 





Em 



:i2) 



where H L = H R is the operator on the LHS of eqn(|5|). 

With the above Hamiltonian there is no mixing between L and R states and hence 
no C = ± splitting. To obtain such a splitting we need a non-zero probability for a L 
state to turn into an R state and vice-versa, which clearly requires some off-diagonal 
terms to appear in the Hamiltonian H. So we alter eqn(|l2|) to define our 'direct mixing' 
Hamiltonian as 



' H L 


a 




a 


Hr 





(13) 
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where we shall choose to keep a real and constant. A simple motivation for such a 
mixing arises from the observation that when the radius of the flux loop is less than the 
flux tube width, we have something that is no longer a distinct loop, but is rather some 
kind of "ball" which will no longer have any definite orientation. In a path integral 
picture, we can think of a path where a loop of orientation L (for example) shrinks into 
a ball, at which point it loses any memory of its initial orientation, and then expands 
back out into a loop of either orientation with equal probablility. This transition will 
lead to a finite amplitude between L and R loops. 

We shall return later to ask how well this model fits the spectrum. For now we 
concentrate on its qualitative predictions, taking the approximation Hl = Hr = M 
where M is the mass of the unmixed state. Then the energy eigenstates are clearly 

M c= ± = M±a (14) 

so that C = ± states are split equally from their common Isgur-Paton value. Three 
immediate observations follow from the above: 

• For states of approximately equal p, the splitting should be roughly the same inde- 
pendent of the phonon number. This would be the case, for example, for the lightest 
J = 0, 1,2 states. However, as we have seen in eqn(^), the splitting does in fact vary 
a great deal amongst these states. 

• It is the average of the C = + and C = — masses that equals the mass with no 
mixing. As we decrease N from iV = oo it is this average that should extrapolate to 
the SU(iV = 2) value of the + mass, because the Hamiltonian there is the same as 
Hl or Hr. (All this up to 0(1/N 2 ) corrections.) However we have seen in Section [11 



that this is not the case: the SU(2) + mass equals the SU(iV > 3) ++ mass up to 
0(1/N 2 ) corrections. 

• We expect the wavefunction to have a smooth limit as N — > oo, and so the probability 
for p to be less than the flux tube radius should also have a non-zero limit. Thus the 
L <r-> R mixing and the consequent C — ± splitting should be non-zero at N = oo, 
just as is observed. (This appears to contradict the conventional statement that 
mixings vanish as N — ► oo but we believe that the standard arguments do not apply 
to our kind of 'mixing'.) 

We might suppose that the first observation becomes irrelevant when we perform the 
actual numerical calculations, but unfortunately it remains an issue. Moreover, however 
we tune the parameters in this model, it remains the case that it is the average of the 
++ and masses that is predicted to continue smoothly to the SU(2) + mass, and, 
as we have seen, this is contradicted by the lattice calculations. Thus other, perhaps 
less straightforward mechanisms need to be considered. 

4.3 Adjoint string mixing 



We pointed out in Section [0]that in SU(iV > 4) theories there exist extra stable strings 



and hence extra states. Since these are just scaled up versions of our fundamental 
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string spectrum they would not help in splitting C = + from C — — . However there 
is another type of string in the theory that we have not yet considered: the one that 
carries adjoint flux. Such a string carries no arrow: it is intrinsically C — +. So any 
mixing would affect the C = + spectrum: the ++ mass would (probably) be driven 
down while the would be left undisturbed. But this does not lead to problems 
with the + in SU(2) because the same adjoint string exists in SU(2) and would drive 
down the + mass there as well. So we have a mechanism that might at last explain 
why the SU(2) + smoothly interpolates onto the ++ for SU(iV > 2). 

The problem with this mechanism is, of course, that the adjoint string can be broken 
by gluon pair production so it is not clear if it makes sense to use a closed adjoint loop 
as the basis for a set of states. However we know that as iV — > oo the string becomes 
stable [[y| so at least for large N it can be used in this way. Since, as we have 
seen, lattice studies || find that the low- lying SU(2) spectrum differs only by small 
corrections from SU(iV = oo), it is reasonable to assume that the adjoint string has 
developed only a modest decay width in SU(2). (As indeed seems to be the case [^31 in 
3+1 dimensions.) If so then the decay time will be long compared to the characteristic 
time scale of the low-lying string modes, and we can safely quantise the string in the 
Isgur-Paton fashion. We will assume that this is so from now on, although an explicit 
lattice verification would clearly be very helpful. 

We can use an extension of the formalism in the previous subsection to derive the 
Hamiltonian, defining the wave-function as 



However since ip a always has C 



H, 




(15) 



adj 








H 







H n 



(16) 



is quite simple: we can clearly reduce it to a two-component calculation in the C = + 
sector, and a simple Isgur-Paton calculation in the C = — sector. We shall assume 
that there is no mixing between fundamental and adjoint loop states that have differing 
phonon occupation numbers. 

In eqnflTGj) H + = H is the usual Isgur-Paton Hamiltonian. H a will be identical 
except that the scale is set by the adjoint string tension, a a , rather than by the funda- 
mental a. It is frequently speculated that a a and o are related by the ratio of quadratic 
Casimirs 

(17) 



a 



N 2 -l 



Lattice calculations in D = 2 + 1 [p2] find that for SU(2) a a ~ 2.5a, which is quite 



close to the value of 8/3 one obtains from eqn([L7|). Thus we expect states based on the 
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adjoint loop to be about a factor of V2.5 ~ 1.5 heavier than corresponding states based 
on the fundamental loop, and so quite massive. Indeed, if we assume a simple two body 
mixing, with H + replaced by M = m — and H a by M a ~ 1.5M, and if we choose a 
so as to obtain the observed m ++ mass, we find that the lightest scalar 'adjoint' state 
has M a Q ++ ~ 1.8m — . This is light enough to be important by its explicit presence in 
the spectrum, in addition to affecting the fundamental states through mixing. 

How can the fundamental flux loop mix with an adjoint loop? Once again our 
underlying physical picture is that, for small p, these loops of string become rather like 
'balls' of flux instead, which may reemerge as a loop of different orientation, or indeed 
even of a different type of flux. We now have the following qualitative predictions: 



As with the direct mixing of Section [4.2|, we expect the mixing, and hence the 



splitting, to be leading order in N, as is observed. 

• As remarked above, since the adjoint loop mixes with just the C = + sector and 
does so for all SU(iV > 2) gauge theories, we expect the SU(2) + mass to continue 
smoothly on to the SU(iV > 2) ++ masses, as is observed. 

• Whether the splitting decreases with increasing phonon number M at fixed radial 
number ur is however not clear - it requires a detailed calculation. 

The observed pattern of the C — ± splitting has proved very constraining on possible 
dynamics, but the model in which the fundamental loop mixes with an adjoint loop 
appears to have the right qualitative features. In the next section we shall see how well 
it can do at the quantitative level. 



5 Fitting the lattice spectrum 

We have performed a large variety of comparisons with the lattice spectrum. We have 
used the different C — ± splitting mechanisms described in this paper; we have fitted 
to all the lattice masses or just to a more reliable subset as discussed above; we have 
used the actual errors in the \ 2 or expanded errors as described below; we have used 
explicit short distance fudge factors of various kinds or a p-dependent string tension as 
described earlier, or indeed no suppression at all. The reader will be relieved to learn 
that we do not intend to describe this very large number of model fits here but will 
rather focus on two of the most relevant examples. Some different model fits can be 
found in || and others will appear in |L0] . 



The generalised flux tube model typically has two or three unknown parameters. 
First there is 7, the (renormalised) curvature/elasticity. Secondly there is the mixing 
parameter characterised by a strength ao- There is the parameter / which characterises 
the short- distance cut-off imposed in the Isgur-Paton model. Alternatively we can use 
a e ff(p) instead of a at short distances and this enables us to do without the parameter 
/. We can then solve the model for various values of these parameters and find which 
fits the lattice spectrum the best. 

In this section we will describe what we find when we fit the lattice data with 
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adjoint string mixing, as described in Section [4.3| , or with direct mixing, as described 
in Section £l|. In both cases we shall use a p-dependent string tension, as described 



in Section 3.3, to embody the short distance corrections to the flux tube picture, in a 



way that introduces no new parameters. In our physical picture of the mixing, we see 



it as occuring at small p and so we choose the mixing term in eqns(|T3|, jig ) to satisfy 



a(p) = a (l-^H 

V a eff 



oo 



(where we now drop the subscript on a a ). In the case of adjoint mixing we determine 
the adjoint string tension from eqn(p!7|). Thus there are two parameters to be fitted: 



7 = 1e — 13/12, where 7# is the string curvature described in Section p72| , and a, the 
mixing strength. 

There are however some problems with just taking the mass spectrum as given in 
Table |l] and doing a least-x 2 fit on the model parameters. First we have already seen 
that the model prediction is going to be far from the lattice 0~ + , so including it in 
the fit might badly distort the final 'best fit'. Moreover there are theoretical reasons 
for thinking that the lattice h spin assignment might be mistaken and that this is 
actually a 4~ + . (And, similarly, that the lattice +_ is really a 4 H .) To a lesser extent 
similar questions arise with the purported J = 1 states. We shall therefore exclude 
these states from the fit but will instead quote the values for these masses, as predicted 
by the best fit to the other states. 

A second problem is that the most accurate masses are for the lightest states. Having 
the smallest errors these will provide the most important component of the x 2 function 
that determines the best fit. On the other hand the lightest states possess the smallest 
radius p and so we expect the flux loop model to suffer the largest corrections for such 
states. Thus we might not want them to dominate, and so possibly distort, the best 
fit. 

This second problem has no unambiguous resolution. The fact is that we know that 
the model is a simple approximation that can at best incorporate only the essential 
features of glueball structure. The most that we can expect is that it roughly reproduces 
the spectrum and whether it can do so is what we want to learn when we fit the model 
to the data. To do this we need to embody what we mean by 'roughly' in some specific 
way into the fitting procedure. In order not to be unduly biased by the very small 
errors on the lightest masses, we have enhanced all the errors by 5% of the mass, and 
add it in quadrature with the statistical error. In practice the best fits turn out to be 
similar whether we perform such an enhancement or not. That is to say, we are saying 
that we want to know if the model can fit the masses within ±5%. Of course the best 
fit might do better than that; but at least it will not be driven by the very small errors 
on the lightest masses for which it is probably least reliable. 

In each case we obtain the predicted spectrum by solving the coupled set of differ- 



ential equations, using the numerical technique detailed in [10|, on a grid in parameter 
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space, and then we perform a final simplex minimization starting at the most favorable 
grid point to find the best x 2 fit- (With just two parameters such a crude approach 
works reasonably well.) The best fit is determined using the lightest three states in the 
++ and sectors, and the two lightest states in the 2 ±+ and 2 ±_ sectors. 

We have remarked that there are reasons to think that these simple models should 
work best in the iV — > oo limit, and so we display the results of fitting to the lattice 
spectrum as extrapolated to iV — > oo. A more thorough display of the results for all 
N may be found in [10]. In Table |] we show the mass spectra for the adjoint mixing 



mechanism, and in Table |3| we present the values of the fitted parameters, for all values 
of N. In Figure |^ we display the spectrum at N = oo. The corresponding spectra for 
the direct mixing mechanism are shown in Table || and in Table |5] are the parameters 
for the direct mechanism. Figure |4| displays the direct mixing spectrum at N = oo. 

We see from the Figures that the best fits, whether obtained using the direct or the 
adjoint mixing mechanisms, are of an overall reasonable quality. A closer examination 
of the detailed spectra does however show that both models have some minor difficul- 
ties. The direct mechanism has problems with the excitations in the J=0 sector: by 
the time the third excitation is reached, which has a greater average p than the ground 
state, the conjugate mixing drives the two model states much farther apart than they 
ought to be. The way that the adjoint mixing model avoids this problem is that the 
mass of the lightest adjoint loop state is naturally close to that of the +_K , because 



Ca/cf ~ 1-5- Thus, after mixing, the ++ * can be largely an adjoint loop, and the 
q++** can then be (largely) the first fundamental loop excitation. However, the direct 
mixing mechanism does have numerically better values of x 2 , and thus we ought to 
consider it more likely than the adjoint mixing mechanism. This suggests that a flux 
tube model that combines both direct and adjoint mixing should be able to do much 
better than either model alone. Combining these models is quite natural; our picture 
for the way an oriented fundamental loop may evolve into an adjoint loop (through 
contracting into a small 'disoriented' ball) is precisely the way we saw the direct mixing 
between fundamental loops of opposite orientation proceeding. This picture will, in 
general, require two mixing parameters. While it is interesting to explore these ideas 
iTif the analysis would clearly benefit from improved lattice calculations where a larger 



number of excited states are accurately determined. 

Returning to the best fit parameters listed in Tables [3] and [5], we observe that all 
our fits require 7 to be positive or very close to zero. Ignoring the values for SU(2), 
there is a trend in the values of 7 consistent with a relationship; see |T(| for more 
details. Taken together with eqn(|6|) this tells us that the observed mass spectra do 
indeed require a non-zero curvature (elasticity) term, 7^ e [0.5,1.0], in the effective 
string model for the confining flux tube. 

From either Figure |] or Figure |] we see that the conventional spin assignments of 
the heavier states are called into question. The model consistently puts the 4~ + and 
the 4 H states at masses corresponding to lattice states with spin 0. These states do 
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not contribute to the x 2 > an d so their masses are purely predicitve-the agreement is 
remarkable. Furthermore, while the lattice states 1 ±+ agree with the model's predic- 
tions for the lightest state with M = 5 (which is the smallest value of M giving spin 1), 
the states conventionally labelled 1 ±_ are at the mass predicted for states with M = 3. 
These states are also related by a spin ambiguity of moduloA on a cubic lattice, and so 
we must also investigate these assignments in the future as well. 

Finally we remark that the features we have described are not only robust against 
the detailed fitting procedures used but much the same conclusions are obtained if we 
replace the p dependence of the string tension with a short distance modification of 
the kind shown in eqn(|j). 



6 Conclusions 

In this paper we set out to test the idea that glueballs are quantised closed strings 
of colour-electric flux. Such a picture arises naturally in linearly confining theories, 
such as SU(iV) gauge theories in 2+1 and 3+1 dimensions, where distant fundamental 
charges are connected by flux tubes. 

We started out with the specific dynamical framework of the Isgur-Paton flux tube 
model |J, in which the excitations of the closed flux loop are either radial or phonon- 
like, and we confronted its mass spectrum with the rather accurate mass spectra avail- 
able in D=2+l SU(iV) gauge theories 0. As N — > oo the gauge theory simplifies in 
ways which bring it closer to some of the model's assumptions, e.g. the neglect of de- 
cays, and so one can argue that a comparison in this limit makes particular sense. If we 
express the observed masses in units of the observed string tension, then the model's 
predictions for these dimensionless ratios involve no free parameters at all. We found 
that these predictions were, for the most part, quite remarkably good in the C = + 
sector of states; and for the C = — sector they embodied the main qualitative features 
even if the quantitative comparison was less good. 

Of course, when the flux tube radius, p, is smaller than the flux tube width the 
picture must break down, embodied in the model by suppressing the potential energy 
below p ~ 1//, where / is a parameter that needs to be determined but which we 
expect to be 0(y/a). In this paper we described other ways of including such a cut- 
off; in particular through a dependence of the string tension on p. By taking er(p) 
from calculations in the literature [[nj, one can avoid having the additional parameter 
/ to fit. We then pointed out that one should in general include a string curvature 
term, which for a closed loop will make a contribution 7_e/p to the effective potential 
that is of the same form as the Casimir string energy. This introduces a parameter 

7 that needs to be determined. Finally we noted that the C = ± degeneracy in the 
model mass spectrum is contradicted by the splittings seen in the lattice spectrum 
and we were compelled to consider dynamical mechanisms that might reproduce these 
splittings. Such mechanisms typically involve a mixing parameter a that also needs 



13 



to be determined. In searching for the best fit to the observed spectrum, we found 
that whatever mixing mechanism we used we invariably required a substantial positive 
string curvature contribution, j E G [0.5, 1]. 

The qualitative features of the observed C — ± splittings proved to be very constrain- 
ing. The only mechanism that we were able to construct that was natural, simple and 
had the right qualitative behaviour, involved adding to the basic flux tube model a 
sector of states built on closed loops of adjoint flux. These are intrinsically C = + and 
we introduced a mixing between these states and the C = + states built on the fun- 
damental loop. We pictured the mixing as arising at small p where a closed flux tube 
becomes a flux-less 'ball', and we conjectured that this kind of mixing may be leading 
order in 1/N, as required by the lattice spectrum. The implication of these calculations 
and also calculations with other mixing mechanisms, such as a direct mixing between 
fundamental loops of opposite orientation, was that the simultaneous presence of ad- 
joint loop and direct mixing was likely to be much more successful in quantitatively 
reproducing the observed spectrum. 

However even these best fits always left us with one very large discrepancy: the h 
state. In the model this is a highly excited state (involving eight phonon units) and is 
predicted to be much heavier than the lattice h . This is a robust result of the model: 
the splitting of the ++ and 2 ++ states, which is two phonon units, is what essentially 
determines the _+ mass. On the other hand the predicted 4 h mass is very close to 
the lattice 0~ + mass (and also for the model's 4 + ~ and lattice + ~). This suggests that 
the lattice calculation may have mislabelled this state; after all one cannot distinguish 
J = from J = 4 by the rotational symmetries of a square lattice. This has provoked 
lattice calculations that are attempting to resolve this rather basic question [T3L |j~5|j. 



(For a detailed exposition of the problem, as well as a means for its solution, see |10|.) 

If there are states built on the adjoint loop they will be about half as heavy again as 
the corresponding states built on the fundamental loop and so only the very lightest 
are likely to be present in the currently available lattice spectrum. Moreover, since the 
adjoint loop can break, these states may have large decay widths, their masses may be 
shifted from their naive values, and perhaps only the lightest states will actually exist. 
However, with a modest improvement in the quality of the lattice calculations, one 
could search for their presence. The same improvement would allow us to search for 
degeneracies between states of differing J but with the same number of phonon 'units'. 
For example there will be an excited ++ with phonon content n+ =2 = n~ l=2 = 1 that 
should be degenerate with a 4 ++ with phonon content n+=2 = 2, 7V =2 = 0. Such (near) 
degeneracies provide a characterisitic pattern that would test the general dynamics of 
the Isgur-Paton flux tube model. 

In SU(iV > 4) gauge theories there are additional strings than the fundamental which 
will be stable if their string tensions are low enough. In D=3+l this is expected to 
be the case on theoretical grounds |16[, and indeed is known to be the case for SU(4) 



17 , |24| . Thus the string model predicts that, if we neglect mixings and decays, the 
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observed mass spectrum will contain towers of states that are exact scaled-up replicas 
of the spectrum arising from the fundamental string, and that the number of these 
'towers' of states will grow as N — > oo. This is a dramatic and robust prediction of the 
basic flux tube picture which can and needs to be investigated by lattice calculations. 

We have seen that the kind of exercise undertaken in this paper, testing a model 
against lattice calculations, has a fruitful impact in both directions. We have been 
forced to generalise the model in ways that, in retrospect, are entirely natural. And 
the model points to both potential weaknesses in the lattice calculations and motivates 
specific further calculations that promise to be very informative however they turn out. 
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Appendix 



In this Appendix we briefly describe how we derive the canonical variables for the 
Isgur-Paton model. First we discuss the traditional model's equation, and then we 
describe how we generalise the formalism to accommodate an effective string tension 
that varies with the loop radius. 

1. The Isgur-Paton Hamiltonian 

We have a circular loop of mass p = 2irp which moves in an effective potential 
provided by the phonon modes etc. as discussed earlier in this paper. Thus the kinetic 
and potential terms are: 

p1 M + t 

T=-p- ; V = 2Kap + F{p)^-^. (19) 
Aixop p 



where F(p) is the usual short distance fudge-factor, which we shall now set to unity 

dp 



for convenience. Under the transcription p p — > we have (with K — 1) 



2 -. 8 8 

T = — — > — dp9p 

2p Ana p ' 1 ' 

We want canonical variables £, P^ such that T is quadratic in In passing from the 

above classical system to the quantum system there is an ambiguity how to order the 

a_ 

various terms in T. We choose P p = -^f so that T oc P p . That is to say we choose 



so that 



t,=P>, (21) 



T^ZL(1) 2 ^ = ^L^ (22) 
Aixo \2J d£ 2 167rcr<9£ 2 v ; 



and we arrive at eqn@. 
2. Generalising to o(p) 

If the vacuum of non-Abelian gauge theories is a (type-II) dual superconducter, then 
flux tubes will arise through a dual Meissner effect. Thus it is interesting and relevant 
to ask how closed flux tubes behave as the radius is varied in the Ginzburg-Landau 
theory of type-II superconductors. As pointed out in what happens is that the 
effective string tension cr e fr(p) vanishes as p — > 0, roughly as in eqn(0). If we replace 
a in the Isgur-Paton model by such a cr c s{p) we will be providing the model with a 
short- distance cut-off that is both natural and employs no free parameters. The price 
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for this is some complication in the quantisation of the model. However, as we shall 
now see, this is a problem that can be overcome. 

As we saw earlier, the tricky term in quantizing the Isgur-Paton model is the kinetic 
energy, so that is where we will begin. Abusing notation slightly, 



(23) 




where a e s(p) is now included in the p dependence of the kinetic energy. If a e s were a 
constant, we would get T as in eqn (|22|) . Now we must find a canonical variable £ such 
that 



d(d_ 
dp dl; 



d 



(24) 



I 1 

<% = (T&(j>)p*dp (25) 
Substituting the function in eqn(^) for cr c fj, we find 

i = a^ f dpp^ (l - e- 1 - 72 ^ . (26) 



where the integration constant has been determined by the boundary condition £(0) = 
0. This integral cannot be put into closed form, yet we need the inverse function 
explicitly so as to substitute it into the potential V(p(£)). While it would be nice to 
include explicitly the function from [ISfl, any function which approximates it reasonably 
well will suffice. So rather than using eqn(^) we shall use a less natural form, but one 
that will suit our manipulations better. Now, we need a function for cr e ff(p) which goes 
to zero as p — > and which approaches a as p — > oo, and which is integrable with 
respect to the measure p^dp. As an example consider 

cr cS (p) = ° (l " e~ /P§ ) 2 (27) 
which, with a suitable value for /, will crudely approximate the dependence found in 



P9f . Performing the integral, and imposing £(0) = 0, we have 

i = a^ I dpp^ ( 1 - e~f^\ =l (T \(p%-\ + -e-f^ ) . (28) 




We cannot find the functional inverse explicitly, but as we are ultimately working with 
a discrete set related to {i,}, we can simply solve the equation numerically to 
give pj which we can put into the potential V = V(p(£j)). (In practice we have used 
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<7 e ff = <t(1 — e ^ p ) 2 which leads to a slightly more complex relation than the above.) 
Returning to the kinetic energy, 

-1 d 2 

(29) 



4tt 0f 2 ' 



and we proceed with the numerical solution of the discretised equation as described 
earlier. 
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SU(oo) 


++ 


4 718(43) 

t: . 1 X U 1 t:U 1 
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t: . V V.' tj \ '.J '-J J 


Q++* 
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6 52 (9) 


6 38(13") 


6 20(13) 
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Q++** 


8 15(15) 


8 23(17") 


8 05(22) 
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7 99(22) 


(1 




6 48(9") 
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6 03(18") 
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o — * 
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Q ** 
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9.19(29) 
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9.47(116) 


2++ 


7.82(14) 


7.13(12) 
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7.19(20) 


6.88(16) 


2++* 






8.51(20) 


8.59(18) 




2-+ 


7.86(14) 


7.36(11) 


6.86(18) 


7.18(16) 


6.89(21) 


2 -+* 




8.80(20) 


8.75(28) 


8.67(24) 


8.62(38) 


2~ 




8.75(17) 


8.22(32) 


8.24(21) 


7.89(35) 


2 — * 




10.31(27) 


9.91(41) 


9.79(45) 


9.46(66) 


2+- 




8.38(21) 


8.33(25) 


8.02(40) 


8.04(50) 


2 +-* 




10.51(30) 


10.64(60) 


9.97(55) 


9.97(91) 


1++ 


10.42(34) 


10.22(24) 


9.91(36) 


10.26(50) 


9.98(25) 


1-+ 


11.13(42) 


10.19(27) 


10.85(55) 


10.28(34) 


10.06(40) 


\ — 




9.86(23) 


9.50(35) 


9.65(40) 


9.36(60) 


1+- 




10.41(36) 


9.70(45) 


9.93(44) 


9.43(75) 



Table 1: Glueball masses in units of the string tension; in the continuum limit ||. The 
SU(oo) values are obtained by extrapolating the SU(iV < 5) values with an 0{1/N 2 ) 
correction. 
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TP modpl 


0++ 
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3.12 


0++* 


6.18(13) 


6.46 


Q++** 


7.99(22) 


8.72 


2 ±+ 


6.88(16) 


6.79 


2 ±+* 


8.62(38) 


9.06 


o-+ 


9.02(30) 


13.86 


4 ±+ 




9.64 


1 ±+ 


10.00(25) 


10.84 


3±+ 




8.30 


0~ 


5.91(25) 


3.12 


* 


7.63(37) 


6.46 


Q ** 


8.96(65) 


8.72 


2±~ 


7.94(35) 


6.79 


2 ±-* 


9.62(66) 


9.06 


o+- 


9.47(116) 


13.86 


4 ±- 




9.64 


1±- 


9.38(60) 


10.84 


3 ± " 




8.30 



Table 2: Glueball masses in units of the string tension. Predictions of the simple 
no-parameter Isgur-Paton flux tube model compared to the actual spectrum of the 
SU(iV = oo) theory. 
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adjoint loop mixing 


group 


a 


7 


Xl/d.o.f 


SU(2) 


0.86±0.65 


-0.19±0.07 


0.47 


SU(3) 


4.5±0.60 


0.57±0.09 


1.1 


SU(4) 


4.2±0.66 


0.42±0.08 


1.2 


SU(5) 


3.8±0.84 


0.25±0.08 


1.4 


SU(oo) 


3.3±0.8 


0.07±0.07 


0.99 



Table 3: Best fit parameters of the adjoint mixing model. 
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4.09 
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7.03 
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Q++** 


8.05 
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8.88 
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8.51 
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5.83 


5.60 


5.36 


Q — * 
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8.70 


8.53 


8.35 
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10.90 


10.77 


10.63 


10.48 


o-+ 


14.06 


14.13 


14.00 


13.93 


13.78 


o+- 




14.21 


14.07 


13.98 


13.82 


4-+ 


11.05 


10.25 


10.13 


10.17 


10.16 


4+- 




11.93 


11.77 


11.59 


11.41 


2 ±+ 


7.55 


6.87 


6.70 


6.70 


6.64 
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9.96 


9.11 


8.97 


9.01 
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8.59 


8.40 


8.19 


7.97 


2 ±-* 




10.83 


10.68 


10.52 


10.35 


1 ±+ 


9.95 


9.19 


9.05 


9.07 


9.05 
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10.87 


10.70 
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10.33 


3 ±+ 


8.78 


8.06 


7.91 


7.92 


7.88 


3 ± " 




9.76 


9.58 


9.39 


9.18 



Table 4: Best fit spectrum of the adjoint-mixing flux tube model to the SU(iV) glueball 
masses, in units of the string tension. 
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group 


a 
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X5/ d - - f 
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-0.42±0.12 
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SU(3) 
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Table 5: Best fit parameters of the direct mixing model. 
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Table 6: Best fit spectrum of the direct-mixing flux tube model to the SU(iV) glueball 
masses, in units of the string tension. 
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Figure 1: Some of the observed SU(iV) C — ± splittings plotted versus 1/N 2 : the mass 
difference between the and the ++ (•) and that between the 2 ±_ and the 2 ±+ (o). 
As N — > oo the dependence is expected to be linear in l/N 2 , i.e. like the straight lines 
added to the plot to guide the eye. 
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Figure 2: The SU(2) 0+ mass (*), the SU(iV > 3) 0++ masses (•), the SU(iV > 3) 0~ 
masses (o), and the average of the ++ and masses (o), plotted against 1/iV 2 ; with 
the expected large- N linear dependence shown in each case. 
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Figure 3: The spectrum for the adjoint mixing mechanism compared to lattice data at 
N = oo. The x-axis gives the total phonon number M. The 0~ + and + ~ are compared 
with M=4, and states with J=l are compared with M=5. 
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Figure 4: The spectrum for the direct mixing mechanism compared to lattice data at 
N = oo. The x-axis gives the total phonon number M. The 0~ + and + ~ are compared 
with M=4, and states with J=l are compared with M=5. 
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